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Abstract 

In a previous paper, we have already considered the system composed by a two 
level atom interacting with a coherent external electromagnètic field. No applica- 
tion whatsoever has been made of the rotating wave approximation. Being specially 
interested in the problem of higher harmònic generations for the case of very intense 
làser fields, we have developed in this letter a much more emcient way to obtain 
these solutions as well as to carry out some calculations in a range in which the pa- 
rameters take extreme vàlues. Also the formalism allows us now to provide analytic 
expressions in the WKB regime for the elèctric dipole moment and the population 
inversion. The spectrum can be decomposed in periòdic and non-periodic contribu- 
tions. Only the latter depends upon the Floquet exponent and can be responsible 
of the main complexities of the observed Rabi revivals and the hyper-Raman shift. 

PACS Numbers: 03.65. Ca, 03.65.Fd and 42.50.Hz 
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In a previous paper we have explored the appearance of higher harmonies as an effect 
arising from the monochromatic elèctric field of very intense làser interacting with a two 
level atom. An adequate theoretical description of the interaction of the atom with the 
strong external field allow us, using the results presented in [Q], to build up the time 
evolution operator of the system and exactly calculate the instantaneous atòmic dipole 
moment and the population inversion. We obtain the Fourier transform of these quantities 
and identifying the present freqüències and its amplitudes we can explain the composition 
of the spectrum of the outgoing radiation. In the present letter an entirely new viewpoint 
will be adopted sheding some light not only on the theoretical aspeets already discussed 
but also enabling us to compute the correspondent physical quantities -either exactly or 
approximate- in a way much more accessible to experimental verification. In order to 
follow a self-contained approach we shall be presenting first a very brief account of (T| . 

Let us consider the physical system containing a two-level atom and its elèctric 
dipole interaction with the coherent field of a làser ||T|. Let \a > and \b > be the atòmic 
eigenstates with energies are and Tiuq respectively. The external monochromatic field 
frequeney is tu and fl is the Rabi frequeney corresponding to an interaction energy hQ . 
In terms of the operators: 

J Q = ^{\b >< b\ - \a >< a\} J_ = \a >< b\ J+ = \b >< a\ (1) 

the system can be described by the well known Hamiltonian 
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H (t) = Uuj ü 



J » + 2 



+ Mïo[J+ + J-]cosut (2) 



Here H(t) is a hermitian element of the triparametric SU (2) Lie- Àlgebra and the physical 
states of the atom can be obtained by acting on the the eigenstates of J with the time 
evolution operator (See refs. : 

U (t) = exp{-^í} exp{r](t)J + } exp{ 7 (í) J } exp{- V *(t) J_} ex P {ih (t) J } (3) 
where 7(t) = Ln(l+ | r](t) | 2 ) and: 

h (t) = -uj t + 2Ü / cos (us) Re [rj (s)] ds (4) 

Jo 

= u* q (t) + 2iq(t) 
u) q(t) - 2iq(t) 

The operator U(t) satisfies the Schròdinger equation if the complex function q(t) is a 
solution of the following second order ordinary linear differential equation (ODE). 

2 

q{t)-2in cos{cüt)q{t) + ( ^-q{t) = ; q (0) = 1 ; g(0) = ^ (6) 
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The main result that we would like to point out here is that it is in fact possible to 
condense the full dynamics of the physical problem in this linear differential equation [p]] . 
Its solution determines all the mathematical quantities that will be considered of interest 
for the physical discussion of the laser-atom interaction. For instance, the rescaled Dipole 
Moment D(t) and the Population Inversion W(t) are defined from the complex solution 
q(t) = r(t) exp{i(fi(t)} of (6) as: 

D(t) = \ q (t)\ 2 -l = r 2 (t)-l W(t) = --r\t)4>{t) (7) 

Using (6) and (7) we can easily check that the following system of differential equations 
holds for D{t) and W{t): 

D{t) + ulD{t) = 2iü íl cos(üüt)W(t) ; D(0) = ; D(0) = (8) 
on 

W{t) = -—co8(ut)D(t) ; W(0) = -1 (9) 

and after some manipulations the following fundamental invariant can easily be con- 
structed: 

D 2 (t) + u 2 D 2 (t) + Lü 2 W 2 (t) = üü 2 (10) 

which guarantees a bound solution for any set of parameters (lüq ^ 0). The differential 
system (8)-(9) has already been discussed in the literature (See Refs. and p). We 

have shown in QXj the relationship of this differential system with the ODE (6), which 
has been proven to be crucial for establishing the non trivial physical features to be 
discussed below. In particular two important properties of this differential system must 
be emphasized: The existence of the first integral (10) and the direct relationship of W(t) 
and D(t) with a complex function q(t) = r(t) exp{i(p(t)} which turns out to be a solution 
of an ordinary linear differential equation with periòdic coefficients. 

The description of the physical system in terms of q(t) instead of the dipole moment 
D(t) and the population inversion W(t) is extremely advantageous as it allows us to 
express directly physical quantities such as the atòmic radiation emitted by the excited 
atom, its spectral composition, the phases as well as the amplitudes of each component 
and the correspondent relationship with the relevant parameters of the system (i.e. the 
atòmic frequency transition u)q, the làser frequency uj and the Rabi frequency VLq). At 
this point we should mention the pioneer papers by Shirley ||, Zel'dovich 0, Cohen- 
Tannoudji and Haroche || and Eberly and coworkers [H] which can be considered the 
first serious attempts to solve this problem rigurously. There are however many aspects 
which have been either overlooked or treated in a different manner by these authors. The 
results hereby presented allows us to analize in a much more efhcient way other interesting 
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effects. Among these we would like to mention the effect of the initial atòmic state as well 
as that of the initial phase of the làser field with or without modulation terms. In order 
to achieve all these goals it seems essential to us a throughout analysis of the properties 
of q(t) and hence those of D(t) y W(t) [[§]. To this end we shall construct the Fourier 
spectrum of the relevant physical quantities. We can also present several approximate 
formulae for the phase and the amplitude for all present modes of the atòmic emission 
spectra. Rescaling the time variable in the form x = ujt and using primes to denote 
derivatives with respect to x, the linear equation (6) reads: 

e 2 e 
q"(x) - 2i-fCOSxq'(x) + —q(x) = q(0) = 1 ç'(0) = i- (11) 

where the two new dimensionless parameters 7 = — and e = — are the ratio between 

r 1 UI UI 

the interaction and foton làser and the atòmic and foton làser energies respectively. 

The ODE (11) can be solved by applying Frobenius Theory pj. Let us define two 
independent functions u(x) and v(x) satisfying u(0) = l,u'(0) = 0, v(0) = 0,v'(0) = 1. 
The Taylor coefficients are given in One can show that: 

• The general solution of (11) is a linear superposition of u(x) and v(x). In particular, 
we set henceforth q(x) = u(x) + i^v(x) 

• The equation posseses a first integral \q'\ 2 + ^-\q\ 2 = C Q which allows us to set 
some bounds for the solutions u(x) and v(x) and yields as a consequence of the first 
integral (11) the absolute bound: 

\u 2 (x)\ + j\v 2 (x)\ = l (12) 

• The set of functions u(x) and v(x) also satisfy: 

u'(x) = --e 2 ^ sinx v*(x) v'(x) = e 2i ^ sinx u*(x) (13) 

These relationships clearly show that u(x) is 27r-periodic (27r-antiperiodic) if and 
only if v(x) exhibits also the same properties. Thus, for any set of parameters 
(e 7^ 0) for which a solution of (11) exists with the property of being 27r-periodic (2ir- 
antiperiodic) we can claim that all solutions will also be 27r-periodic (27r-antiperiodic). 
Notice also that u(x) and v(x) have the following symmetry properties: 

e 2 

u(x + 7r) = u(tï)u*(x) — —v*{tí)v*{x) (14) 



v(x + 7r) = v(ir)u*(x) + u*(tt)v*(x) 



(15) 
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Also u(— x) = u*(x) and v(—x) = —v*(x). Clearly these properties can be used 
to know the vàlues of the functions u(x) and v(x) for all x if they are known just 
on the [0,7r] interval. This interval may be reduced in fact to [0, -]. In particular, 
§(u(2tt) + v'(2n)) = Re[u(2n)] = 1 - 2(efle[«(f K(f )]) 2 . 

• The set of functions D(x) and W(x) may be expresed as: 

D(x) = elm(u(x)v*(x)) W(x) = -Re[ e - 2i ~< smx (u 2 (x) + -v 2 (x))} (16) 

• As the equation (11) has periòdic coefficients a Floquet analysis can be performed. 
According to this theorem any solution can be expressed in the form: 

q{x) = Ae ivx F + (x) + Be~ ivx F_ (x) (17) 

where F+(x) y F_(x) are 27r-periodic functions. Using the above properties we end 
up with an explicit expression for the Floquet exponent, namely: 

v = -L arccos[ite[u(27r)]] = - arcsin[ei2e[|M(~)i>*(£)|]] < v < i (18) 

Z7T 7T 2 2 2 

Its physical significance lies on the fact that it actually determines the spectral 
decomposition of the emitted atòmic radiation. The two extreme cases correspond 
to solutions 27r-periodic (u = 0) or 27r-antiperiodic iy = |). 

The Fourier spectrum of u(x) and v(x) takes the following genèric form: 

j=oo j=oo 

f{x) = AJ2 Fje i{3+u)x + B F^e i( j-"> (19) 

j = — OD j = — OQ 

where the Floquet exponent v and the Fourier coefficients Fj can easily be found 
simultaneously with the help of the following recurrence relation: 

^{(p + v) 2 -j}F p = ( p + u + l)F p+1 + {p + u-l)F p ^ (20) 

which can be solved through a numerical method using continued fractions Q. The 

instantaneous dipole moment of the atom can finally be expressed as: 

D(x) = elm{u{x)v*{x)) = A x (x) + A 2 (x) (21) 

This expression is in fact the sum of two different contributions. The A\(x) term is a 
superposition of odd harmonies of the làser frequeney. The A 2 (x) term is a superposition 
of even harmonies shifted up and down an amount given by 2u. 

oo 

A 1 {x) = ^D iC os[(2i + l)x] (22) 

3=0 

oo oo 

A 2 (x) = ^Dfcoa[2{j + v)x]+Y i Djcoa[2{j-v)x] (23) 

j=0 3=1 
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The correspondent Fourier amplitudes are: 

M 2 - M 2 ~ ~ 
^ = 2 (Me 2 + M J )2 E F r F r+23+l M e = F 2r (24) 



r=— oo r=— oo 

oo oo 



^ = - 2 (M f;M oT r E (-1)^-,^ M °= E F W (25) 



r=— oo 



The Fourier spectrum of D(x) dictates the composition of the emitted atòmic radiation in 
interaction with a làser. A qualitative picture shows the appearance of a triplet centered 
in the odd harmonies with freqüències (2s + 1 ± 5)cu. The amount of the shift 5 is given 
by the Floquet exponent (18) in the form 5 = 1 — 2u and coincides with the generalized 
Rabi frequeney. The non-periodic contribution A 2 (x) of the instantaneous atòmic dipole 
moment, which can be interpreted as a superposition of odd harmonies shifted up and 
down an amount 5, is the origin of the hyper-Raman peaks of the triplet. Likewise, the 
population inversion can also be cast in the form: 

W(x) = -1 - 2- r coszD\z)dz = ILÍx) + U 2 (x) (26) 
e Jo 

as the sum of two contributions Hi(x) (a superposition of even harmonies of the làser 
frequeney) and H 2 (x) which is a superposition of odd harmonies shifted up and down an 
amount given by 2u: 

oo 

U^x) = ^Wj coa[2jx] (27) 

3=0 

oo oo 

n 2 (x) = J2 w ^ cos [( 2 j + 1 + 2u ) x ] + J2 w f cos [( 2 j + 1 - 2u ) x ] ( 28 ) 

j=l j=0 

with Fourier amplitudes given by: 

oo 

W = -(1 + W ) - + W; + W. ) (29) 

3=1 

Wj = -?(£),- + + Dj ~ £>J ·- 1 ) ; j>\ (30) 
W - = -](DÏ + Dï + °L·^£) (31) 
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Aside from the numerical procedure provided by (18), (24)-(25) and (29)-(32), a com- 
plete analytical reconstruction of the spectrum can be carried out starting with an input 
provided by the D(x) and W(x) functions. If one can find analytic expressions for the 
functions u(x) and v(x), the vàlues for D(x) and W(x) are inmediately known from (16) 
and the amplitudes Dj, Df, Wj, and Wf can easily be found as: 

4 fi 

Dj = -J A x (x) cos[(2j + l)x]dx (33) 
2 f~ 

Df = - / 2 A 2 (x)cos[(2(j ±u)x]dx± 

J 71 JO 

1 — 

± . - . / 2 {A 2 (x-tt)- A 2 (x + 7T)}sm[{2{j±u)x]dx (34) 

Wj = Tli(x) cos[2 jx}dx (35) 

2 f~ 

Wf = - I U 2 (x) cos[(2j + 1 ± 2v)x\dx T 

J 71 JO 

=F ——— [ 2 {Tl 2 {x-n) -n 2 (x + 7r)}sm[(2j + l±2v)x}dx (36) 
ix sm(27rz/) Jo 

Where m =1 and rrij =2 (J > 1). These expressions are obviously exact as long as so 
were the knowledge of the Floquet exponent v and the functions u(x) and v(x). 

As a concrete example of the above discussion, we shall be considering the case in 
which the parameters of the physical system verify 7 w e >> 1. This condition holds 
for vàlues of the làser frequency smaller than the vàlues of the resonance frequency and 
the coupling constant at = - is of the order of one. In this way we are considering 
subresonant systems uj « u ~ Q with interaction and transition energies of similar 
vàlues. Since in other papers 0|,[§J a numerical study of this range of parameters has 
been made we have chosen the same range for the sake of comparison. This case seems 
also quite adequate for using the WKB method for high vàlues of the parameter e. In 
regard to the equations (11)-(13) this means that we are looking for a solution of in the 
form q(x) = exp{i^ sin x + e^^lo z k(x)e~ k }, where the Zk(x) verify: 

z' 2 (x) = ~(1 + 4a 2 cos 2 x) (37) 

4(z) = -^^{4-1^) + £ 4(^)4-^)} *>i (38) 

ZZ 0V X ) j=l 



After calculating the sequence for higher orders of approximation, u(x) and v(x) can be 
expressed as: 

u{x) « V ^ { (l--)A(a,6,x) + (l + -)/ 2 (a,6,x)} (39) 

2 1 + 4a 2 cos 2 x)4 A A 
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v(x) ~ ^{f2(a,e,x)- fl (a,e,x)} (4Q) 
eVA (1 + 4a 2 cos 2 x) 4 

fi(a,e,x) ps [(A + 2a) (\A + 4a 2 cos 2 x - 2a cos x)]2$(a, e, x) (41) 

f 2 (a,e,x) ps [(A - 2a) (a/Ï + 4a 2 cos 2 x + 2a cos x)] ^ $* (a, e, x) (42) 

r/ eAMx,À; 2 ] (l + 8a 2 )Mx,A; 2 ] -Fix,*; 2 ],, 
$(a,e,x) ps expW ^ + " LL· - 1 )} (43) 

where A = a/1 + 4a 2 , A; 2 = -^"^ and F[x, A; 2 ] and E[x, k 2 ] are the incomplete elliptic 
integrals of the iirst and second kind respectively. iFrom these expressions one can derive 
the vàlues of the Floquet exponent as well as the vàlues of the instantaneous dipole 
moment. These vàlues take the form: 

«fca) = I(^l + ' 1+8 ° 2 »ffl- A '^) (44) 
7r 6eA 

v(e, a) ps arcsin[| siníl(e, a)|] < z/ < - (45) 

For sufficiently large vàlues of e and practically any value of a (e >> a) the Floquet 
exponent v is a 2e G - periòdic function of e of the form: 

v ~ 7T- */ < e < e (46) 

i/ w e <e<2e (47) 

where e = 7r(2>/ï + 4a 2 £'[ 1 ^ 2 ])~ 1 and f repeats itself every time that e increases in 
2e . For a given value of the coupling these oscillations do not die away and keep going 
even if the làser frequency decreases. In particular when a = 1 (the energy of interaction 
equals the energy of the transition) the value of eo is approximately 0.596086. Notice 
the predictive power of (44) and (45). One can use it for selecting vàlues of e and a with 
a given Floquet exponent v. Thus the spectral composition of the dipole moment and 
the population inversion can be selected at will with a judicious choice of the parameters. 
Finally one can easily conclude that: 



{x , L 2E[k 2 ]((l + 8a 2 )E[k 2 ]~K[k 2 ]) ] (4g) 



2 VI + 4a 2 £[À; 2 ] V 3r 2 7r 2 

If r ps n, v = and the functions u(x) and v(x) are 27r-periodic. Both the dipole moment 
D(x) and the population inversion W(x) are composed by pure even and odd harmonies 
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of the làser. If r ~ n + |, v — |, and -u(x) and are 27r-antiperiodic. The spectrum of 
-D(x) contains just odd harmonies and that of W(x) just even harmonies. If r ~ + |), 
v — \i and the functions -u(x) and are 47r-antiperiodic. The atòmic dipole moment 
can be constructed in terms of its following two components: 

. . . 2a cosx , . 

Ai(a,x) = -- - — (49) 

A vl + 4cr cos z x 



2a cos{e\E[x, fc 2 ] + ^{(1 + Sa^Ejx, fc 2 ] - F[x, fe 2 ]}} 

A 2 (a,e, x) = — (50) 

A v 1 + 4a 2 cos 2 x 

Notice that Ai(a,x) is independent of e and it is a purely even 27r-periodic term under 
the periòdic translation Therefore, this term enjoys the adequated symmetry to represent 
the periòdic part of the dipole oscillation. One can easily identify a term like this as an 
oscillation at the làser frequeney, with modulated amplitude and just composed by odd 
harmonies of this frequeney. The term A 2 (a,e,x) has quite different properties being 
periòdic just for those vàlues of the parameters for which the Floquet exponent v takes 
rational vàlues. The vàlues of the instantaneous dipole moment hereby found can be used 
to find its respective spectral components. The periòdic part yields the odd harmonies 
spectrum with amplitude and phase are given by: 

8a rleosxcos^ + l).] 

3 vrVl + 4a 2 Jo Vl + 4a 2 cos 2 x V ' 

For sufficiently lower làser freqüències (sufficiently large e) the odd harmònic spectrum 
results quite independent of this frequeney. The different amplitudes depend just on the 
strenght of the laser-atom interaction and smoothly decrease for higher order harmonies. 
In particular the amplitude and phase of the first harmònic can explicitely be given by: 

^-J^lï^l-ï^ï^ (52) 

The second contribution gives rise to the hyper-Raman spectrum. For any set of the 
physical parameters of the system, either íl + v or Q — v must be an integer number. One 
can can construct the correspondent amplitudes in the form: 

D+ = 4a t\ cosfcA + ±g£)E[x fc 2 ] - ± K [x, fc 2 ] ± 2{j + u)x] ^ Q 
3 7rA Jo y/1 + 4a 2 cos 2 x 



Aa ri cos[( £ A + ±&-)E[x, fc 2 ] - ±K[x, fc 2 ] T 2(j - u)x] j>1 
7rA Jo y/1 + 4a 2 cos 2 x 
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where the upper sign holds for + v an integer and the lower sign for the integer being 
Q — v. A large number of oscillations substantially lowers the value of the integral, 
the correspondent vàlues of the amplitudes leads one to consider easily the r egion of 
the spectrum of maximal amplitudes. Therefore for the first case, the amplitude of the 
harmonies with freqüències 2j + 2v will be negligible while the harmonies with freqüències 
2j — 2v will have significant amplitudes which can even have similar intensities to that of 
the odd harmonies. These large amplitudes will be distributed over a region with vàlues of 
2j close toQ + v. This behaviour justifies the typical harmònic generation spectrum found 
by experiment showing an intensity distribution that first decreases rather steeply as the 
harmònic order increases (the odd spectrum) and then after remaining almost flat for a 
number of harmonies finally decreases steeply again in a "plateau" form (the non-integral 
spectrum) |ID[ . All these facts are in fairly good agreement with the cases discussed in 
H where the shift 5 is given by an elliptic integral as in our case, but the correspondent 
intensities are given in terms of Bessel functions. Likewise, the population inversion can 
explicitely be characterized by using the two contributions: 

n x (a,x) = - 1 === (55) 

Ay 1 + 4cr cos z x 



Aa 2 cosx cos[(eA + ±g£)E[x, k 2 } - ±F[x, k 2 }} 

n 2 (a, e, x) = , — — (56) 

A v 1 + 4a 2 cos 2 x 

and as in the previous case Ui(a,x) is independent of e. It is a purely 7r-periodic term 
with the adequate symmetry to represent the periòdic part of the population inversion 
whose Fourier spectrum is just composed by even harmonies of the làser frequeney with 
amplitudes: 

2rrij /"§ cos[2jx] 



W j = . 3 „ / ; 1 J 1 = (57) 

3 tt\/1 + 4a 2 Jo y/1 + 4a 2 cos 2 x 

Where m =1 and rrij =2 (j > 1). The term n 2 (a, e, x) depends on the two characteristic 
parameters of the system and yields the rest of the spectrum with amplitudes: 

8a 2 fí cos x cos[(eA + ^f)E[x, k 2 } - ^K[x, k 2 } ±(2j + l + 2v)x\ 



W+ = / : — b^j. — 1 bex i ' — j — ; /_i -j>o 

3 7rA Jo \/T+ 4a 2 cos 2 x 



n- = / f cos x cos[(eA + ^)^h fc2 ] ~ é; K i x > k2 ] I gj + 1 - 2v)x\ . v i 

7rA Jo vT+ 4a 2 cos 2 x ' _ 

with the same convention of signs as the one given above. 
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Conclusions. In a previous paper the present authors had already considered the 
system composed by a two level atom interacting with a coherent external electromagnètic 
field by solving the main equations in the Schròdinger picture without using the rotating 
wave approximation. We were specially interested in the spectral composition of the 
atòmic dipole moment and the population inversion. These quantities are determined by 
means of two functions u(x) and v(x) which are independent solutions of an ordinary 
differential equation with periòdic coefficients. Floquet analysis applied to this equation 
yields the freqüències which are present in the spectra and allows us to identify two 
different contributions in D(x) (the atòmic dipole moment) and W(x) (the population 
inversion) which yield the different spectral components through quadratures. In this 
letter we have presented analytical expressions for the freqüències and amplitudes 
within a particular range of vàlues of the physical parameters. These expressions 
may be used to analyze the influence of the different parameters in the spectra. Our results 
are also clearly of experimental interest in assigning odd and Floquet-shifted character to 
each observed harmònic. This gives rise to the decreasing and flat parts of the observed 
spectrum. 
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